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ABSTRACT 


Beamforming  of  the  TFO-37  array  reduces  the  rms  of  the  noise 
up  to  14  db  over  an  average  single  sensor,  and  the  signal/noise 
improvement  approaches  14  db  also  depending  on  the  band  pass  fil¬ 
ter  used.  The  signal  loss  after  beamforming  and  summing  is  approx¬ 
imately  1  db..  Power  spectra  was  also  performed  on  the  noise  after 
summation  as  well  as  on  individual  traces  and  shewed  a  15  db  re¬ 
duction  at  1  cps ;  this  is  an  improvement  of  N  over  an  average  sin¬ 
gle  element. 

Maps,  digital  programs,  instrument  frequency  response,  etc. 
are  included 
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INTRODUCTION  AND  PURPOSE 


The  purpose  of  this  report  is  to  present  a  preliminary, 
evaluation  of  the  new  TFO  array.  This  array  consists  of  37 
JM  short  period  vertical  seismometers  in  the  form  of  three 
concentric  hexagons  of  six,  twelve,  and  eighteen  sensors 
each  with  one  in  the  center.  The  approximate  diameter  is 
thirty  kilometers. 

In  addition  to  presenting  the  results  of  a  basic  inves¬ 
tigation,  another  intention  is  to  compile  all  information 
pertaining  to  the  array  into  one  report  which  could  be  used 
as  a  reference  for  future  studies. 

With  this  in  mind,  a  number  of  appendices  are  included 
consisting  of  the  following: 


Appendix  A: 

Map  of  Array 

B: 

Seismometer  description  and  locations 
in  latitude  and  longitude,  etc. 

C: 

Digital  tape  format 

D: 

Digital  programs  written  for  the 
evaluation 

E: 

Calibration  parameters 

F: 

Sample  output  from  which  Table  1  was 
prepared 

G: 

JM  frequency  response 

The  primary  objective  is  to  determine  the  amount  of  sig¬ 
nal  loss,  noise  attenuation,  and  S/N  gain  which  result  from 
beamforming  digital  recordings  of  four  teleseismic  events. 

The  LASA  Bulletin  statistics  for  these  earthquakes  are  shown 
in  Table  2.  Since  this  report  includes  such  a  small  sampling 
of  events,  one  from  each  quadrant  was  selected  to  give  full 
azimuthal  coverage. 
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PROCEDURE 


Digital  seismograms  are  recorded  at  the  site  on  7-track 
magnetic  tape  in  multiplex  form  at  the  rate  of  556  bpi.  Sine 
wave  calibrations  (1  cps)  are  performed  daily  on  each  data 
channel . 

Each  seismogram  used  in  this  study  is  detrended,  thereby 
removing  mean,  and  each  was  converted  from  digital  counts  to 
equivalent  earth  motion  (mu)  at  1  cps.  For  comparative  pur¬ 
poses  processing  was  accomplished  twice  using  different  band 
pass  filters.  One  has  a  cut-off  at  .4-3  cps,  the  common 
LASA  or  SDL  filter,  while  the  other  cuts  off  at  .7  -  5  cps 
and  will,  hereafter,  be  referred  to  by  the  designation  Fc. 

Both  fall  off  at  the  rate  ox  about  24  db/octave  (Figure  1) ; 
however,  since  the  phase  angle  of  the  data  was  preserved, 
the  effective  decay  rate  becomes  48  db/octave.  As  can  be 
immediately  seen,  the  Fc  filter  removes  considerably  more  of 
the  low  frequency  microseisms,  the  effect  of  which  will  be 
perceived  in  the  results. 

The  seismograms  were  then  time-shifted  to  align  P  arrivals 
by  applying  computed  travel  time  differences.  No  attempt  was 
made  to  fit  by  any  other  method. 

Selected  portions  of  the  multiplexed  data  amounting  to 
two  minutes  of  noise  and  one  minute  of  signal  were  demultiplexed 
and  converted  to  SDL  library  format.  The  resulting  data  were 
combined  under  one  label  and  input  to  the  beamforming  program. 
Calibrations  were  also  converted  to  the  format  under  several 
labels  but  were  fed  directly  to  the  magnification  program. 
Outputs  from  this  program  were  saved  on  punched  cards  and 
printouts.  Quality  control  has  been  limited  primarily  to  the 
elimination  of  dead  f>r  reversed  traces. 
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RESULTS 


Table  1  contains  the  final  results  of  the  amplitude  anal¬ 
ysis  and  contains  three  major  headings  under  Signal,  Noise  and 
S/N.  Although  the  table  is  self-explanatory,  a  few  of  the 
salient  features  can  be  mentioned.  Relative  to  the  signal  there 
is  a  slight  loss  in  amplitude  between  the  mean  and  the  sum  as 
expected  possibly  due  to  small  mis— alignments i  In  addition,  the 
apparent  db  loss  can  also  be  attributed  to  the  signal/noise  ratio. 
The  Tunisia  earthquake  is  poorly  recorded  and  the  "p"  probably 
contains  some  noise  which  precludes  an  accurate  estimate  of 
signal  amplitude  during  the  phasing  process.  A  deterioration  in 
amplitude  (mu)  is  obvious  between  the  SDL  and  Fc  filter  by  about 
4  db  in  all  cases.  This  suggests  that  the  SDL  filter  is  more 
effective  in  preserving  signal  amplitude  within  the  frequency 
range  of  the  events  studied.  On  the  other  hand,  an  examination 
)f  Figure  2  reveals  the  undersirable  loss  of  first  motion  as 
well  as  a  loss  of  signal  character  caused  by  the  SDL  filter. 

The  SDL  filter  appears  to  over-simplify  the  signature  at  the 
same  time.  This  retrogression  is  due  to  the  low  frequency  of 
the  recorded  signal  which  is  well  down  in  db  from  the  .4  cps 
low  cut-off. 

With  respect  to  the  noise,  both  filters  have  an  acceptable 
reduction  in  rms  but  the  Fc  is  the  better  of  the  two.  This  is 
due  to  the  removal  of  the  low  frequency  microseisms  which  pre¬ 
dominate.  The  amplitude  of  the  maximum  noise  excrusion,  in¬ 
dicated  by  the  heading  RANGE,  is  also  shown.  The  ratio  of  the 
Range/rms  is  approximately  3/1,  a  value  normally  considered 
when  referring  to  peak  vs  average  rms . 

The  last  heading  in  the  table  is  probably  the  most  crucial, 
statistically  speaking,  since  it  deals  with  the  all-important 
S/N  ratio.  Again,  both  filters  function  well  but  the  Fc  responds 


best  by  about  3  db  over  the  SDL .  However,  neither  equals,  on 
the  average,  the  /~N  criterion. 

Again,  regarding  the  Signal  portion  of  the  table,  it  may 
be  found  that  alignment  of  the  P  phase  is  best  achieved  either 
by  using  average  observed  travel-time  anomalies,  or  by  phasing 
using  cross-correlations. 

Power  spectra  was  also  performed  on  the  noise  of  the  four 
events  by  a  digital  program  based  on  the  Blackman  and  Tukey 
method.  Figures  3  and  7  illustrate  the  results  of  spectral 
analysis  on  100  seconds  just  prior  to  the  signal  of  the  phased 
sum  of  all  channels  vs.  the  average  spectra  of  outputs  recorded 
by  18  to  19  spatially  selected  individual  elements.  It  was 
felt  that  taking  the  average  of  a  number  of  sensors  was  some¬ 
what  better  than  selecting  just  one  and  that  spatial  sampling 
of  the  array  would  give  about  the  same  results  as  processing  the 
total. 

As  to  be  expected,  the  spectral  shape  is  approximately 
the  same  for  any  given  filter  and  will  not  vary  too  much 
whether  the  analysis  is  performed  on  an  individual  trace  or 
the  phased  sum  of  many.  In  all  examples  shown,  summing  the 
elements  attenuates  the  noise  at,  for  example,  1  cps  approx¬ 
imately  15  db  over  the  average  single  instrument  which  is 
equal  to  a  factor  of  N.  The  best  S/N  ratio,  again  referring  to 
Tabla  1,  is  several  db  lower.  For  comparative  purposes.  Figure 
3  shows  the  results  of  raw  data  processing  of  the  Hokkaido 
event.  This  reveals  that  the  lower  frequencies  are  basically 
coherent  across  the  array  and  are  preserved  even  after  time 
shifting.  Where  normally  the  peak  power  occurs  at  just  under 
1  cps  after  filtering,  the  peak  on  raw  data  is  at  about  0.25 
cps . 
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SUMMARY 


The  results  of  this  basic  study  revealed  several  prominent 
features  : 

1)  The  P  phase  from  strongly  recorded  events  is  attenuated 
only  slightly  by  beamforming.  However,  it  appears  that 
the  weaker  the  event  the  greater  the  loss. 

2)  The  rms  of  the  noise  is  reduced  by  11-14  db,  depending 
on  the  filter,  after  summing. 

3)  The  zero-to-peak  maximum  noise  excursion  in  mu  is  about 
3  times  that  of  the  average  rms  obtained  in  a  100 
second  window. 

4)  There  is  also  an  11-14  db  improvement  in  S/N  ratio  by 
beamforming . 

5)  There  is  an  approximate  15  db  noise  reduction  at  1  cps 
as  shown  by  the  power  spectra. 

6)  Neither  filter  used  is  completely  satisfactory.  The 
SDL  seems  to  preserve  signal  amplitude  at  the  expense 

} 

of  character  but  does  not  reduce  the  noise  to  a  min¬ 
imum.  The  Fc  filter  reduces  the  noise  almost  to  the 
system  background,  but  attenuates  the  signal  amplitude 
while  preserving  the  first  motion  and  shape. 


RECOMMENDATION 


As  a  result  of  this  study,  three  areas  of  interest  were 
found  which  should  be  investigated;  they  are: 

1)  The  use  of  average  observed  travel-time  anomalies 
and  cross-correlations  to  reduce  beamforming  sig¬ 
nal  loss. 

2)  The  design  and  evaluation  of  a  set  of  band— pass 
filters  to  determine  the  most  efficient  with 
respect  to  signal  preservation  and  noise  reduction. 

3)  The  comparison  of  TFO-37  with  LASA,  NORSAR  and  the 
other  arrays  in  terms  of  their  respective  detection 
and  identification  capabilities. 
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RELATIVE  MAGNIFICATION 
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Figure  3.  Power  spectra  of  unfiltered  noise 
Hokkaido,  25  February  1968. 
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Figure  4.  Noise  power  spectra,  Hokkaido, 
25  February  1968. 
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Figure  5.  Noise  power  spectra,  South  Pacific, 
2A  February  1968. 
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Figure  7..  Noise  power  spectra,  Tunisia, 
25  February  1968. 
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APPENDIX  B 


Seismometer  descriptions  and  locations 
Tonto  Forest  Seismological  Observatory 


TOWTO  rOREST  SEISMOLOGICAL  OBSERVATORY  -  37  ELEMENT  ARRAY 
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APPENDIX  C 


Digital  tape  format  * 


The  TFSO  digital  recordings  are  written  on  7-track  IBM  compatible 
magnetic  tape. 

1.  DEFINITIONS 
1 .  1  DATA  WORD 

A  36-bit  word  containing  three  samples  of  data  from  three  different 
seismographs.  Sources  (seismometers)  are  identified  as  channel  "i" 
where  i  is  the  detector  identity  and  1  £  i  s  48. 

1 . 2  TIME  WORD 

A  36-bit  word  which  contains  clock  samples  in  milliseconds  and  the  time 
the  clock  was  reset  to  zero  from  0000Z. 

1 . 3  SCAN 

Equal  to  sixteen  data  words,  i.e. ,  one  data  sample  from  each  of  48 
seismographs. 

SCAN  =(16  data  words)  (3  samples  per  data  word)  =  48  samples. 

The  real  time  required  for  a  scan  is  approximately  50  milliseconds.  The 
scans  alternate  starting  times  to  synchronize  the  input  with  the  output.  The 
scans  are  started  at  intervals  of  48  and  52  milliseconds  in  the  following 
sequence: 

0,  48,  100,  148,  200,  248,  300,  348,  400,  etc. 
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1.4  DATA  RECORD 


Consists  of  one  Time  Word  and  220  Scans.  All  the  data  for  all  channels 
for  a  11  second  sample  period  (220  data  samples  from  each  seismometer). 
Data  Record  =  (220  scans)  (16  data  words  per  scan)  +  (  1  time  word)  = 

3,  521  thirty-six  bit  words. 

1 . 5  HEADER  RECORD 

An  auxiliary  record  containing  only  36  bits.  It  is  the  first  record  of  each 
tape  (If  there  is  more  than  one  file,  it  is  the  first  record  of  each  file.) 
The  data  on  this  record  includes  the  date,  format  designation,  and  the 
A strodata  amplifier  gain  settings. 

2.  FORMATS 

2.1  OVERALL  FORMAT 


»4 

v  c, 
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V  V 

SC  Pi 


T3 
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w 


A  Tape  can  record  approximately  2  hours  and  20  minutes  of  data  with  each 
data  record  being  11  seconds  of  data. 


lTA  record 


Record  Gap 


Time  Word 


Data  Wordi 


Data  Word- 


Scan  1 


Scan  2 


Data  Word^ 

ch  46 
ch  47 
ch  48 

ch  1 

Data  WordjY 

ch  2 

ch  3 

ch  4 

Data  Wordjg 

ch  5 

_ chJa. 

Data  Word. 


ch  46 

ch  47 
ch  48 


Scans  3 
thru  219 


Start  Time  of  first  scan 


Data  Word 


?  505 


Data  Word 


3506 


Scan  220 


ch  46 

Data  Word,_,.ch  47 
_ ch  48 

Record  Gap 

Next  Record 


FORMAT  OF  ONE  DATA  WORD  CONTAINING  THREE  DATA  SAMPLES  FROM  THREE  DIFFERENT 


2.4  FORMAT  OF  HEADER  RECORDER: 


01234  5  6  78910111213  14  1516171819  20  21  22  23  24  25  26  27  28  29  30  31  32  33  34  35 


loooooo 

0000  0  0  0  0  0 

00000  0|  000000 

oooooooool 

^ - : - . - -j 

.  A  J 

V. _  J 

Format  Astrodata  Spares  Year  Day  Year 

884  Amplifier 
gain  settings 


2.5  FORMAT  OF  TIME  WORD: 


0  1  2  34  56  78  9  10  11  12  13  14  15  16  17  18  19  20  21  2223  2425  26  2728  29  30  31  32  33  34  35 
0  0  0000000  OjO  Q  Q  Q  Q  Q  Q  Q  Q  Q  OOOQoooooonnnnn 

sv.  A: 


5 


u 

(T) 

a 

tn 


B 


A  -  Number  of  five  minute  B  -  Number  of  ms  from  the  time  in  A  to  the  first 

intervals  past  0000Z  when  data  sample  of  the  record.  The  value  of  the  low 

the  clock  was  reset  to  zero.  order  bit  is  21  instead  of  the  standard  2°. 

Therefore,  the  bit  value  of  B  is  2^,  2 23  2^ 

23,  22,  2 1 .  . 

The  decoded  time  word  gives  the  time  of  the  start  of  the  first  scan  of  each 
record . 


CHANNEL  ASSIGNMENT 


Seismometer 
Channel  Assignment 


Seismometer 
Channel  Assignment 
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♦Station  Timing  System 


APPENDIX  D 

Digital  programs  written  for  the  evaluation 

SEISMIC  DATA  LABORATORY 
Alexandria,  Virginia 
Digital  Computing  ‘Section 

A.  INDENTIFICATION 

Title:  Mergence  of  2  Seismograms  Under  One  Label 

COOP  Identification:  Z148  MERGSEIS 

Category:  General 

Programmer:  P.A.  Gillespie 

Date:  September  30,  1967 

B.  PURPOSE 

I  his  program  will  take  any  two  seismograms  from  one  tape 
and  merge  them  togehter  under  a  new  label  with  50  channel  ident¬ 
ifiers. 

C.  USAGE 

lm  Operational  procedure:  This  is  a  Fortran-63  main  pro¬ 
gram  with  one  subroutine.  The  subroutine  calls  two  routines 
from  the  Monitor  Library  pack.  They  are:  (1)  SKIPREC  and  (2) 
SKIPFILE.  A  data  card  is  needed  for  each  seismogram.  If  the 
seismogram  number  is  the  same  for  both  seismograms,  two  data 
cards  must  still  be  used  with  the  same  number  on  each  card.  The 
input  and  output  tapes  may  be  either  library  or  subset  format. 

2 •  Input  card  format: 

Card  requesting  the  first  seismogram. 

Parameter  Meaning 

Seismogram  # 

Type  of  input  tape 
Type  of  output  tape 


(a) 

Column 


1-10 

11-13 

14-16 


LABEL 
IN  * 
OUT* 


*The  numerical  value  of  IN  and  OUT  has  no  meaning,  but 
must  be  non- zero .  However,  the  sign  of  this  numerical  value 
determines  the  type  of  tape.  A  positive  number  for  either 
IN  or  OUT  or  both  designates  a  library  tape,  whereas  a  negative 
number  refers  to  a  subset  tape. 

(b)  Card  requesting  the  second  seismogram. 


This  card  is  identical  to 
the  seismogram  number,  which  may  or 

3.  Print-outs : 

(a)  Error  messages. 

Comment 

TYPE  OF  INPUT  TAPE  NOT  SPECIFIED- 
JOB  TERMINATED 

TYPE  OF  OUTPUT  TAPE  NOT  SPECIFIED- 
JOB  TERMINATED 

INTERNAL  ERROR  IN  CALCULATION 
THE  ERROR  MAY  BE  FOUND  IN  AN  IF 
STATEMENT  WITH  THE  NUMBER  (/vV) 

IN  ONE  OF  THE  THREE  BRANCHES 


the  first  card,  except  for 
may  not  be  different. 


Cause 

Cols.  11-13  on  either  one 
or  both  input  cards  was 
left  blank  or  a  zero  used. 

Cols.  14-16  on  either  one  or 
both  input  cards  was  left 
blank  or  a  zero  used. 

A  (996)  in  the  parenthesis 
means  that  either  the  first 
or  second  seismogram  has 
more  than  25  channel  identi¬ 
fiers.  A  (997)  in  the  paren¬ 
thesis  means  that  the  input 
label  does  not  have  the  cor¬ 
rect  number  of  channels  on 
it. 


(b)  Printed  -  Output 

After  the  two  seismograms  have  been  merged,  the  program 
will  print  out  a  message  saying  this  along  with  the  new  label. 
The  new  seismogram  number  and  the  number  of  channels  is  a  sum¬ 
mation  of  the  two  input  seismogram  numbers.  The  channel 


identifiers  are  consectively  in  order,  taking  up  two  output  lines. 
The  fourth  line  of  the  label  is  the  event  identifier,  which  is 
taken  from  the  second  seismogram,  with  the  new  seismogram  number. 

4.  Space  required;  18000 

5.  Caution  to  Users: 

(a)  The  system  must  have  Monitor  Library  Routines 
named  (1)  SKIPREC  and  (2)  SKIPFILE. 

(b)  The  seismogram  numbers  on  the  two  input  data 
cards  must  be  in  the  same  sequence  as  they  are 
on  the  input  tape . 

(c)  If  the  seismogram  number  is  the  same  for  both 
of  the  input  seismograms ,  then  the  same  number 
must  appear  on  both  of  the  input  data  cards. 

(d)  The  sign  of  IN  should  be  the  same  on  both  input 
data  cards. 

(e)  The  sign  of  OUT  should  be  the  same  on  both  input 
data  cards. 

(f)  The  signs  of  IN  and  OUT  need  not  be  the  same. 


SEISMIC  DATA  LABORATORY 
Alexandria,  Virginia 
Digital  Computing  Section 


A.  IDENTIFICATION 

Title:  Data  Converter  for  TFO  Acquisition  System  Tapes 
COOP  Identification:  M224  TFXARY  II 
Category:  Information  processing,  conversion 
Programmer:  J.G.  Ferry 
Date:  2  January  1968 

B.  PURPOSE 

This  program  makes  a  standard  library  seismogram  tape  from 
digital  tapes  in  the  format  of  TFO  acquisition  system  tapes. 

C.  USAGE 

!•  Operational  Procedure:  This  is  a  FORTRAN  and  CODAP 1 
main  program.  A  disk  file  has  to  be  on-line.  If  more  than  one 
input  tape  is  used,  a  pause,  222,  will  occur  with  the  next  tape 
reel  number  on  the  typewriter. 

When  a  pause  occurs  with  zero  in  the  A-register  the 
tape  reel  number  from  output  unit  1  should  be  entered  into  the 
A-register  in  octal. 

2.  Parameters : 

NC,  number  of  channels  to  be  decoded 

ST(1),  hour  of  start  time 

ST (2),  minute  of  start  time 

ST  (3),  second  of  start  time 

SR,  sampling  rate  in  points  per  second 


TL,  length  of  record  in  seconds  to  be  decoded 
IL ( 1 ) ,  seismogram  number 
NREL,  input  tape  reel  number 

(NCNS (I)  , 1=1 , NC)  channel  numbers  to  be  decoded 

(IL  (I) ,1=8, 37,  second  and  third  seismogram  label  cards 

3.  Space  required:  19,571  locations  plus  first  four 
disk  file. 

4.  Temporary  storage:  none 

5 .  Printouts : 

A  message  is  printed  on  the  standard  output  unit 
if  an  end -of-file  is  encountered  on  logical  unit  3. 
Program  continues,  but  uses  length  of  data  available 
instead  of  TL. 

6.  Error  returns:  nor. a 

7 .  Error  stops;  none 

8 .  Input  and  output  tape  mountings ; 

Input  tape  3,  output  tape  2,  output  tape  1,  out¬ 
put  tape  7  (BCD),  output  tape  11  (BCD),  output  tape  24 
(operator  comment) . 

9 .  Input  and  output  formats ; 

A.  Input 

1 .  Cards : 

First  card 

Column 

1-5 
6-10 
11-15 
16-20 
21-30 
31-40 


Format  (I5,3F5.2,F10.3,2I10) 

Data 

NC 

ST  (1) 

ST  (2) 

ST  (3) 

SR 

TL 


41-50  IL (1) 

51-50  NREL 


Second  card 


Format  (251?) 


Column 


Data 


1-3 

4-6 

7-9 


NCNS(l) 
NCNS  (2) 
NCNS  13) 


Third  card  Format  (25A3) 

Data:  second  card  of  seismogram  label,  IL(8)  to  XL (32) 
Fourth  card 

Data:  third  card  of  seismogram  label,  IL(33)  to  IL(42) 

NOTE:  The  four  card  data  deck  is  repeated  for  each 

seismogram. 

2 .  Tape: 

The  input  tape  is  a  digital  tape  from  the  TFO 
acquisition  system  with  a  one  word  header  record 
followed  by  2641  words  per  record  and  one  end-of  file 
at  the  end  of  the  tape.  Each  2641  word  record  must 
contain  a  timing  word  and  220  scans  of  48  channels. 

B.  Output 

1.  Output  tape  on  unit  1  is  a  library  seismogram 
tape. 

2.  Output  tape  in  unit  2  is  a  plot  tape  with  internal 
control. 

3.  Output  tape  on  unit  7  (BCD)  is  a  listing  of  the 
events  on  the  output  tape  from  unit  1. 


4.  Output  tape  in  unit  11  (BCD)  is  the  first  card 
of  the  seismogram  label  and  index  card  for  each 
event  on  the  output  tape  from  unit  1. 

10.  Selective  stops: 

At  pause  222  the  typewriter  will  display  the  next 
reel  number  to  be  mounted  on  input  unit  3,  or  pause  0 
for  entering  the  reel  number  on  output  unit  1  into  the 
A-register  in  octal. 

11.  Time  required; 

Approximately  1/2  minute  per  channel  for  180  seconds 
of  seismogram  length. 

12.  Accuracy : 

Same  as  input. 

13.  Caution  t;o  users; 

Maximum  number  of  channels  per  seismogram  is 
25.  Channel  numbers  (NCNS(I))  to  be  decoded  must  be 
between  1  and  48  inclusive.  Maximum  length  of  seis¬ 
mogram  record  (TL)  is  800  seconds. 

14.  Equipment  configuration: 

l  COOP, 23551, MUNZL, 1/3/  ’/2/l/S/56/57/E/7=51/ll=52/24=54 , 

30,10000. 

15.  References : 

None 

D.  APPENDIX 


Astrodata  Seismic  Data  Acquisition  System  Tape  Format. 


Astrodata  Seismic  Data  Acquisition 
System  Tape  Format 

SECTION  A.  General  Description: 

The  Astrodata  Seismic  Data  is  recorded  on  a  7-track  tape 
with  a  density  of  556  bits/inch.  The  physical  records  on  the 
tape  consist  of: 

Header  record 

IRG  (inter-record  gap) 

Scan  record  1 
IRG 

Scan  record  2 


Scan  record  N, 

EOF  Gap 
EOF 

where  each  scan  record  contains  eleven  seconds  of  multiplexed 
data  for  48  channels  with  a  sampling  rate  of  20  points  per 
channel/per  second.  More  detail  is  given  in  Sections  B  and  C 
below. 

SECTION  B.  Notes: 

1.  Standard  36  bit  word  IBM  format  is  used  with  parity  in 
channel,  row,  number  7.  Parity  is  odd. 

2.  Longitudinal  parity  is  taken  after  the  header  and  each 
data  record. 

3.  Data  samples  begin  at  the  start  of  each  scan  and  consist 
of  two  characters,  12  bits,  each.  Characters  are  numbered  consecu¬ 
tively  in  each  scan:  for  example,  characters  5  and  6  make  up  the 
data  sample  for  channel  number  3  in  each  scan. 


4.  One  data  record  consists  of  21,  126  characters. 

5.  Record  gaps  appear  as  follows: 

Record  gap  (prior  to  record  number  1)  =  0.75  in 
IRG  =  0.75  in. 

I OF  Gap  =  3.5  in. 

6.  EOF  is  standard  -  17 g 


SECTION  C.  Bit  significance  in  Example  Consisting  of — Two_Data 
Records : 

1.  Header  Record  -  one  word-12  octal  number s-written  at 
the  start  of  each  file. 

000/001/000/001/010/000/000/000/101/100/110/011  BINARY 
qj_q120005463  OCTAL 

A  B  C  D  E 

A.  format  designation;  B.  Astrodata  884  amplifier  gain  settings; 
C.  no  significance;  D.  year  designation;  E.  day  of  year. 

2.  Data  Record  -  time  word  written  at  11.000  second  in¬ 
tervals  at  the  start  of  each  data  record.  In  data  record 
number  1  the  time  word  is: 


ooi/ooi/ooo/oio/ooi/oii/ooo/10Q/101/1QQ/OQO/1Q1» BINARY 

'  C 


A 


B 


A.  No  significance 

B.  Three  octal  numbers  designation  the  number  of  five  minute 
intervals  past  0000Z  when  the  clock,  in  C,  was  reset  to 
"zero."  In  this  example,  time  was  reset  at  2208  or  at  1200Z. 

C.  26  bits  designation  the  time  interval  in  ms  from  the  time  in 
B.  to  the  start  of  the  first  record.  Bit  value  is  as  follows: 

226,  225,  224,  223, - 24,23,22,21  ms 

In  data  record  number  2  the  time  word  is  exactly  the  time 
word  in  number  1  plus  11000  ms.  A  scan  is  initiated  every 


50  ms*>  therefore,  the  sample  time  may  be  calculated  as: 

Time  specified  by  the  time  word  +  50  (scan  number  -1)  ms 


Data  -  Data  record  number  1  -  scan  number  1: 

Characters  1  and  2  have  the  following  significance: 
2°  is  equivalent  to  4.9  mv. 


+  2 
1  0 


^^2^2^2^2^2"^2^2^2^2^’2^ 

0001110101 


Channel  number  1  sampled  at  the 


time  specified  in  the  time  word. 

Characters  3  and  4  make  up  channel  number  2  at  the  same  time. 
Similarly  for  all  48  channels  in  first  scan. 

Data  record  number  1-scan  number  2: 

Similar  to  scan  1,  except  the  sampling  time  is  50  ms  later*. 
Data  record  number  1-scan  number  220: 

Sample  time  is  the  time  word  plus  10,950  ms. 


♦Actually,  the  intervals  are:  0 , 48 , 100 , 148 , 200, 248 ,300 , etc. 


SEISMIC  DATA  LABORATORY 

Alexandria,  Virginia 
Digital  Computing  Section 


A.  IDENTIFICATION 

Title:  Analysis  of  TF037  seismological  data 
COOP  Identification:  Z  146  TFOSAN 
Category:  General 
Programmer:  P.A.  Gillespie 
Date:  February  29,  1968 

B.  PURPOSE 

To  process  TF037  data,  20  samples  per  second,  less  than 
4000  point  per  channel  from  a  special  SDL  library  tape.  This 
special  SDL  library  tape  is  made  from  processing  TF037  data 
through  program  M224  TFXARY I I ,  and  then  through  program  Z  148 
MERGSEIS .  All  data  to  be  evaluated  by  TFOSAN  must  follow 
this  procedure.  The  calibrations  to  be  used  as  data  input 
must  first  be  processed  by  Z153  TFOCAL . 

C.  USAGE 

1*  Operational  Procedure:  This  is  Fortran-63  main  program 
with  twelve  subroutines.  The  main  program  calls  five  routines 
from  the  Monitor  Library  Pack.  They  are:  (1)  LIBDATE ,  (2) 

ERASE,  (3)  SKIPFILE,  (4)  DISC  63,  and  (5)  ABMAX . 

2.  Parameters :  The  deck  set-up  for  TFOSAN  to  run  from  a 
binary  deck  is  as  follows: 

(a)  COOP  card  -  3  COOP,  23551,  MUNZL ,  I/1/0/4/6/S/56/57, 

10,1000 


(b)  9  EXECUTE,  4,  50,  1 

(c)  Binary  Deck 

*  (d)  Seismogram  Card  (1) 

*  (e)  Calibration  Cards  (one  for  each  channel  on  the  input 

tape) 

*  (f)  Parameter  Card  (1) 

*  (g)  Filter  Coefficients  cards  (4)  with  the  type  of  filter 

in  cols.  1  -  4  of  the  first  filter  card  on’y. 

*  (h)  Delta  time  and  anomaly  cards  (5) 

*  (i)  Deletion  card  (2)  (purpose  is  to  prevent  a  channel 

from  being  processed.) 

*  (j)  Type  of  input  tape  (1)  (this  card  requires  a  signed 

integer  constant;  if  the  sign  is  positive  it  assumes 
a  library  type  tape,  if  it  is  negative  the  program 
assumes  a  subset  type  tape)  > 

*  See  Section  C-9  "Input  and  Output  formats"  for  the  details 
of  (d)  -  (j),  and  any  pertinent  restrictions. 

2a.  Cover  Card:  Place  in  Special  Operating  Instructions  section 
of  cover  card.  "When  pause  zero  occurs  place,  the  number  of  the 
tape  mounted  on  logical  unit  #6,  in  the  A-register  in  octal 
form. " 

3.  Space  required:  20348 

4.  Temporary  storage  required:  Disc 

5.  Printouts : 

6.  Error  returns ;  Are  self-explanatory 

7.  Error  stops:  None,  except  when  an  error  message  is 
printed  out. 

8 .  Input  and  Output  Tape  Mountings ; 

(a)  Input  Data  Tape  on  logical  unit  #1 

(b)  Output  plot  tape  on  logical  unit  #4 

(c)  Output  save  tape  on  logical  unit  #6  (This  save 
tape  has  all  the  aligned  traces  plus  two  channels, 
the  unphased  and  phased  sums.) 


9.  Input  ar.d  output  formats: 


(a)  Input  tape  (logical  unit  #1) 

The  input  tape  may  be  either  a  standard  SDL  un¬ 
packed  library  tape  with  a  modified  label,  or  a 
standard  SDL  Subset  tape. 

1)  The  modified  label  for  the  standard  SDL  Library 
Tape  has  50  channel  identifiers  instead  of  the 
standard  25  channel  identifiers. 

2)  The  Subset  tape  will  have  the  standard  SDL  subset 
label. 

(b)  Output  tapes 

1)  Output  plot  tape  on  logical  unit  #4 

This  tape  will  plot  on  the  160A  the  aligned 
traces  of  the  raw  input  data  channels,  and  the 
four  sum  traces.  They  are  all  plotted  to  the 
absolute  maximum  in  the  data. 

2)  Output  save  tape  on  logical  unit  #6 

This  save  tape  has  all  of  the  aligned  trace 
plus  only  two  of  the  sum  traces,  the  unphased  and 
phased  sum. 

(c)  Input  data  cards 

1)  SEISMOGRAM  CARD:  (1)  card  Format  is:  (110) 


Column 

Parameter 

Meaning 

1-10 

ISEIS 

To  locate  the  desired 

seismogram  on  a  tape 

2)  Calibration  cards:  There  must  be  one  card  for 
each  channel  on  the  input  tape.  Format:  (7X,F8.0) 


Column 

Parameter 

Meaninq 

1-3 

The  channel  identifier  for 

which  this  calibration  is 

intended 

8-15 

DMAG 

The  calibration 

3)  Parameter  Card:  (1) 

card  Format:  (2A3,  2F3.0, 

F5.1,  V2 

.0,  15) 

Column 

Parameter 

Meaninq 

1-3 

OPTION* 

To  filter  the  individual 

traces  only 

4-6 

OPTION2* 

To  filter  the  sum  traces 

only 

*  The  only  values  that  may  be  used  for  these  two  para¬ 
meters  are  the  alphabetical  words  "YES"  and  "NO". 

7-9  REF HR  The  hour  of  the  earliest 

arriving  signal 

10-12  REFMI  The  minute  of  the  earliest 

arriving  signal 

13-17  REFSE  The  second  of  the  earliest 

arriving  signal 

18-19  XPP  The  number  of  seconds  for 

the  P-P/2  gait.  If  a  zero 
is  used  it  will  automatically 
be  set  for  8  seconds. 


24  INTAPE  The  number  of  the  tape  to  be 

mounted  on  logical  unit  #1 

4)  Filter  coefficient  cards:  (4)  cards  Format  is 

(A4,1X,E15.7,3E20.7/(5X,E15.7,3E20.7) ) 

Column  Parameter  Meanina 


INTAPE 


Column 

1-4 


TYPEFILTER 


6 “20  FILCOA(4  to  a  card) 
21-40 
41-60 
61-80 


Meanina 

The  type  of  filter  being  used 
to  filter  the  data  with. 

FILCOA  are  the  "A"  coefficients 
FILCOB  are  the  "B"  coefficients 


5)  Delta  times  and  anomaly  cards  (5)  cards 
Format  is:  (10(F3.2,  F5.3)) 

Column  Parameter  Meanina 


Column 

Pa  ram 

1-3 

DT 

4-8 

AM 

9-11 

DT 

12-16 

AM 

17-19 

DT 

20-24 

AM 

25-27 

DT 

28-32 

AM 

33-35 

DT 

36-40 

AM 

41-43 

DT 

44-48 

AM 

49-51 

DT 

52-56 

AM 

57-59 

DT 

Delta  time  for  first  channel 
Anomaly  for  first  channel 
Delta  time  for  second  channel 
Anomaly  for  second  channel 
Delta  time  for  third  channel 
Anomaly  for  third  channel 
Delta  time  for  fourth  channel 
Anomaly  for  fourth  channel 
Delta  time  for  fifth  channel 
Anomaly  for  fifth  channel 
Delta  time  for  sixth  channel 
Anomaly  for  sixth  channel 
Delta  time  for  seventh  channel 
Anomaly  for  seventh  channel 
Delta  time  for  eighth  channel 


60-64 

AM 

Anomaly  for  eighth  channel 

65-67 

DT 

Delta  time  for  ninth  channel 

68-72 

AM 

Anomaly  for  ninth  channel 

73-75 

DT 

Delta  time  for  tenth  channel 

76-80 

AM 

Anomaly  for  tenth  channel 

The  next  four 

cards  are  a 

repeat  of  inis  one. 

6)  Delection 

cards:  (2) 

cards  Format:  25  (A3) 

Column 

Parameter 

Meaning 

three  at 

a  time 

I  DUMP 

When  the  channel  identifier  in 

this  field  matches  the  channel 

identifier  in  the  same  place  on 

the  label  that  channel  will  be 

deleted  from  the  processing 

NOTE:  The  channel  identifiers  on  these  two  cards  must  be 
in  the  same  logical  position  as  they  are  in  the  label  in 
order  to  match.  For  example  if  the  first  channel  identi¬ 
fier  in  the  label  is  Zl,  be  sure  to  check  if  the  blank  is 
to  the  left  of  the  Z  or  the  right  of  the  1  and  place  in 
the  same  position  on  the  deletion  card,  as  this  is  an  A3 
format.  Otherwise  they  will  not  match. 

7)  Type  of  input  tape:  (1)  card  Format:  I 10 

Column  Parameter  Meaning 

1-10  INTAPE*  Type  of  input  tape 

*  The  integer  value  in  this  field  has  no  meaning,  but 
must  be  non- zero.  However,  the  sign  of  the  number  de¬ 
termines  the  type  of  tape.  A  positive  integer  designates 
a  library  tape,  whereas  a  negative  number  refers  to  a 
subset  tape. 


D .  METHOD 


(1)  Detrend 

(2)  Demagnify 

(3)  Filter  (optional)  (4-pole  Butterworth  -  16  coeff.) 

(4)  Reverse  filtering 

(5)  Compute  RMS  on  100  secs  of  the  noise  -  (N) 

(6)  Obtain  maximum  peak-to-peak  amplitude  of  the  same  100 
sec  of  the  noise,  and  divide  by  two  N(P-P/2) 

(7)  Obtain  maximum  peak-to-peak  amplitude  of  the  signal 
divide  by  two  (S) 

(8)  Compute  the  ratio  of  the  signal  to  the  RMS  of  the  noise 
S/N 

(9)  Compute  the  ratio  of  the  signal  to  the  per  ■  peak 
divided  by  two  of  the  noise  S/N  (P-P/2) 

(10)  Form  an  unphased  sum  of  each  point  for  all  channels 

(11)  Form  a  phased  sum  of  each  point  for  all  channels  based 
upon  the  earliest  arrival  of  the  signals 

(12)  One  the  unphased  and  phased  sum  channels  are  formed, 
these  channels  are  then  processed  for  steps  3  -  8. 

(13)  Brennan  Weights  are  computed  -J.S//.n2)  ^  .  g  ,  where 

I (S/N) i 
i=l 

NC  =  no.  of  channels;  and  then  applied  to  each  channel. 

(14)  Unphased  sum  trace  of  the  channels  with  the  Brennan 
Weights  applied 

(15)  Phased  sum  trace  of  the  channels  with  the  Brennan  Weights 
applied 

After  the  Brennan  unphased  and  phased  sum  traces  are 
formed  they  are  evaluated  using  steps  3-8. 


(16) 


SEISMIC  DATA  LABORATORY 


Digital  Computing  Section 


IDENTIFICATION 

Title:  Subset  Library  or  Subset  Tapes 

COOP  Identification:  M219  SUBSETSL 
Category:  Information  Processing 

Programmer:  R.L.  Hawkins 

Date:  2  November  1967 


PURPOSE 

To  subset  a  seismogram/or  seismograms  contained  on  the  Four 
standard  SDL  tapes: 

1.  Unpacked  library  tape:  the  127  word  label  record  is 
written  first,  followed  by  the  data  channel  records 
(number  of  channels S-25)  with  up  to  16000  points  per 
channel.  After  each  channel  an  EOF  is  written,  except 
for  the  last  channel  which  has  two  EOF's. 

2.  LASA  library  tape:  the  same  as  the  unpacked  library 
tape  except  the  last  channel  may  have  either  a  double 
EOF  or  an  EOF  followed  by  -1. 

3.  Subset  tape:  contains  a  label  record  of  127  words 
followed  by  the  data  chs-^el  records  and  then  one  EOF. 

4.  Packed  library  tape:  is  written  in  the  same  manner  as 
a  subset  tape  except  the  data  channel  records  are  packed 
as  NC  12  -  bit  values  per  word,  where  l£.NC^4.  That  is, 
for  NC  =  4  we  have, 

12  BITS  12  BITS  12  BITS  12  BITS 

IX  (1)  CHAN  ( 1 )  CHAN  (2)  I  CHAN(3)|  CHAN(4)|pata  pt.l 


IX  (2) 


IX  (N) 


CHAN ( 1 ) 

CHAN  (2) 

CHAN  ( 3 ) 

CHAN (4)  Data  pt.2 

•  • 

•  • 

•  • 

•  • 

•  • 

•  • 

•  • 

•  • 

CHAN (1) 

CHAN  (2) 

CHAN (3) 

CHAN (4) |Data  pt.N 

C.  USAGE 

1.  Operational  procedure:  This  is  a  F-63  main  program  and  a 
CODAP- 1  subroutine  UNPKLTP 

2.  Parameters : 

a.  A  card  containing  the  plot  parameters  -  ICODE  and  RANGE 

b.  A  card  containing  the  seismogram  number,  the  start  time 
of  the  subset,  the  range  or  duration  time  of  the  subset, 

a  plot  switch,  a  tape  switch,  and  a  pause  switch. 

c.  A  card  containing  the  ID-codes,  such  as  SPZ,RH,DW^,etc. , 
of  those  data  channels  desired. 

Cards  b)  and  c)  are  repeated  if  more  than  one  seismogram 
is  requested. 

3.  Space  required:  This  program  is  dimensioned  for  16,340  words, 
a  maximum  of  16000  data  points  /  channel  may  be  requested 
from  any  of  the  four  standard  tapes. 

4.  Temporary  storage  requirements :  none 

5.  Alarms :  none 

6.  Printouts :  For  each  seismogram  request  the  following  is 
printed  out: 

1)  The  output  subset  label 

2)  Scale  factors  for  the  plots  of  all  data  channels.  If  no 
plots  are  requested,  no  scale  factors  are  printed. 

3)  The  absolute  maximum  value  of  each  data  trace. 

7.  Error  printouts; 

1)  If  the  seismogram  requested  is  not  on  the  input  tape,  the 
following  is  printed  out,  "Seismogram  XXXXX  not  on  Tape", 


and  the  next  request  is  processed. 

2)  If  the  same  channel  is  requested  more  than  once,  the 
following  is  printed  out  "channel  XX  is  requested  twice", 
and  the  next  seismogram  request  is  processed. 

3)  If  the  requested  ID-code  does  not  equal  any  of  the  ID- 
codes  (ID  (8) -ID  (32) )  from  tape,  the  following  is  printed 
out,  "channel  XX  not  on  input  tape",  and  the  next  seis¬ 
mogram  request  is  processed. 

8.  Error  stops;  none 

9.  Input  and  output  tape  mountings: 

a)  Input  tape  unit  3 

b)  Output  save  tape  unit  2 

c)  Output  plot  tape  unit  4 

10.  Input  card  formats: 


card  1 

Format 

(FlO ,  110) 

Col. 

Data 

Description 

1-10 

RANGE 

Vertical  height  of  plots  in  .01  inches. 

Be  sure  that  a  range  is  selected  which 
does  not  cause  plots  to  overlap. 

11-20 

ICODE 

Number  of  . 01  inches  between  successive 

data  points  on  plots. 

card  2 

Format 

(110,2 (F4,F3,F5) ,16,110) 

Col. 

Data 

Description 

1-10 

NSELS 

Seismogram  number 

11-14 

BHR 

Beginning  time  of  subset  request  in  hours. 

15-17 

BMN 

Beginning  time  of  subset  request  in  min¬ 

utes  . 

18-22 

BSE 

Beginning  time  of  subset  request  in  sec¬ 
onds  . 

23-26 

RHR 

Range  time  or  duration  of  the  subset  in 

hours . 


27-29 

30-34 

35-40 

41-50 

51-60 

61-70 

card  3 

Col. 

1-3 

4-6 


RMN 

RSE 

LPLOT 


LORS 


IPAUSE 


Format 


Range  time  or  duration  of  the  subset  in 
minutes . 

Range  time  or  duration  of  the  subset  in 
seconds . 

Plot  switch.  If  L  plot  is 

-1-  labeled  plots 
0-  no  plot 

1-  non-labeled  plots 
Subset  switch.  If  LORS  is, 

-1-  subset  of  a  packed 
library  tape 
0-  subset  of  either  an 
unpacked  or  LASA  library 
tape 

1-  subset  of  a  subset  tape 
number  of  channels  to  plot  across. 

A  non- zero  value  causes  the  program 
to  pause  so  another  tape  of  the  same 
type  can  be  mounted.  This  option  is 
used  only  when  requesting  a  subset 
from  tapes  of  the  same  type  (UNPACKED, 
PACKED,  OR  SUBSET) 

(26A3) 


Data  Description 

IDCHAN  (1)  First  channel  ID-code  so,  written  on  the 
label  of  the  input  tape. 

IDCHAN  (2)  2nd  channel  ID-code 


IDCHAN  (NC)  the  channel  ID-code. 

The  following  important  features  are  contained 
within  the  program  primarily  to  enable  the  user  to 
run  SUBSETSL  with  ease.  Some  of  these  features  are 


taken  from  the  SUBSET  program  because  they  were  found  to  be 
extremely  useful  through  usage;  however  a  few  like  the  range 
time  in  hours  were  retained  simply  because  it  was  felt  and 
consequently  discovered  that  many  users  of  SUBSET  had  built 
up  a  conditional  response  that  the  range  time  in  hours  is  in 
columns  23-26  of  the  second  data  card,  even  though  a  range  in 
hours  is  completely  unrealistic  at  any  practical  sampling  rate. 

Important  SUBSETSL  Features 

1)  If  the  beginning  time  field  (columns  11-22  of  card  2)  is  left 
blank,  the  origin  time  of  the  seismogram  on  the  input  is  used. 
For  example  if  the  seismogram  starts  at  10  15  20  Z,  the  sub¬ 
set  will  begin  at  10  15  20  i. 

2)  If  the  range  time  field  (columns  23-34  of  card  2)  is  left  blank, 
the  duration  of  the  subset  runs  to  the  end  of  the  seismogram. 

For  example  if  the  seismogram  ends  at  15  10  00  Z  and  the 
beginning  time  request  of  the  subset  is  15  09  00  Z,  the  subset 
contains  1  minute  of  seismic  data. 

3)  If  the  beginning  time  request  is  less  than  the  origin  time 
of  the  seismogram,  the  subset  will  start  at  the  beginning  of 
the  seismogram.  For  example  if  the  user  requests  a  subset 
start  time  of  05  10  30.0,  but  the  origin  time  of  the  seismogram 
is  05  10  40.0,  the  subset  will  begin  at  05  10  40.0. 

4)  If  the  range  time  is  larger  than  that  which  is  available  on 
the  seismogram,  the  duration  of  the  subset  runs  to  the  end  of 
the  seismogram.  For  example  if  the  seismogram  starts  at  10 

10  10Z  and  ends  at  10  11  10Z  and  at  range  time  of  2  minutes  is 
requested,  the  subset  will  contain  1  minute  of  data. 

5)  There  is  an  automatic  pause  built  into  the  program  to  enable 
the  mounting  and  hence  processing  of  a  different  type  of  tape. 
For  example  if  the  user  is  retrieving  seismograms  from  a 
packed  tape  and  then  wants  to  process  an  unpacked  tape,  the 
program  will  pause  automatically.  However  to  process  two  tapes 
of  the  same  type,  the  pause  flag  (column  70  on  data  card  2) 


must  be  set  to  a  non- zero  value. 


6)  The  tape  type  (pack,  unpack,  or  subset)  is  printed  out  on  the 
typewriter  for  each  seismogram  processed.  This  feature  allows 
the  operator  and  user  to  keep  track  of  what  has  been  processed 
if  program  termination  occurs. 

7)  If  the  user  trys  to  retrieve  rhe  same  data  channel  twice  or 
more,  this  seismogram  request  will  not  be  processed.  Use 
either  the  Z120  RODRDATA  or  Z140  RODBUDSC  to  obtain  the  same 
data  channel  more  than  once. 

8)  Incorrect  specification  of  the  ID-codes,  such  as  punching  DHl 
on  data  card  3  when  you  want  DW1,  results  in  a  waste  of  com¬ 
puter  time  and  termination  of  that  seismogram  request.  Make 

sure  that  what  is  punched  on  data  card  3  is  truly  on  the  input 
tape. 

9)  To  retrieve  every  channel  that  is  written  on  the  input  tape 

tor  a  desired  seismogram,  merely  punch  "ALL"  in  columns  1-3 
of  data  card  3. 

10)  Just  about  any  feasible  kind  of  plotter  spacing  within  the 
limits  of  present  plotting  capabilities  is  allowed, 
a)  If  the  number  of  plots  written  across  is  left  blank 
(column  60  of  data  card  2) ,  two  options  exist; 
i?  Whenever  the  number  of  data  channels  requested  is 

greater  than  10,  ten  channels  are  written  across  on  the 
plot  tape. 

ii;  If  the  number  of  channels  requested  is<il0,  that 
number  of  channels  will  be  written  across  for  plot 
control. 

b.  Specify  the  number  you  want  written  across  on  the 
plot  tape  in  columns  51-60  of  data  card  2. 

c.  If  the  user  is  retrieving,  for  example,  one  channel 
from  each  of  ten  seismograms  and  wants  these  ten  channels 
written  across,  don't  insert  a  10  in  columns  59-60  of 


data  card  2.  Tha  program  control  statement  is  written  on 
tape  unit  4  every  time  a  seismogram  is  processed.  Hence  the 
above  request  results  in  more  plot  expectations  than 
availability  of  plots  to  fill  the  logical  records.  If  the 
user  wants  to  have  the  above  plot  presentation  or  something 
of  the  same  kind  of  logic,  write  plotter  switch  11  on  and 

the  number  of  plots  across  on  the  plot  section  of  the  white 
job  card. 


11)  The  following  label  will  be  written  (punch  -  1  in  columns 
39-40  of  data  card  2)  before  each  set  of  plots  correspond¬ 
ing  to  that  seismogram  request. 

a.  SUBSET  OF 

b.  SEISMOGRAM  NO  XXXXX 

c.  Heading  information  from  seismogram  label  like 
AP-OK  3  JUL  67 

d.  Channel  ID-codes  for  the  plots  like  ZH  RH  TH  T 

12)  The  scale  factors  for  each  plot  is  printed  out. 

13)  The  absolute  maximum  value  of  each  data  channel 
requested  is  printed  out. 

14)  A  127  word  label  for  each  seismogram  is  written  before 
the  data  records.  The  format  is  as  follows: 

IL  (1)  =  seismogram  number 

IL  (2)  =  the  number  of  data  channels 

IL  (3)  =  the  number  of  data  points/channel 

IL  (4)  =  the  sampling  rate 

IL  (5)  =  start  time  of  subset  in  hours 

IL  (6)  =  start  time  of  subset  in  minutes 

IL  (7)  =  start  time  of  subset  in  seconds 

IL  (8)-IL  (32)  =  ID-codes  of  data  channels 

IL  (33)  -  IL(42)  =  Identification  or  pertinent 

information  of  the  seismogram  such  as  station 

code  (i.e.  AP-OK),  date  of  event,  name  of  event, 

and  IL(42)  contains  the  seismogram  number  presently. 


IL (43 ) -IL  {127)  =  not  used  presently 
11.  Selective  jump  and  stop  settings:  none 

1?.  Timing;  Program  running  time  is  dependent  upon  the  following 
factors ; 

1)  The  kind  of  subset  requested 

For  example,  if  a  subset  of  a  packed  tape  is  requested, 
the  fixed  point  data  must  be  unpacked  through  the  use  of 
subroutine  UNPKLTP . 

2)  Where  on  the  input  tape  the  desired  data  is  written 

If  the  desired  seismogram  is  not  at  the  beginning  of  the 
tape,  library  routines  SKIPFILE  and  SKIPREC  must  be  used. 

3)  The  number  of  data  channels  and  number  of  points/channel 

4)  Whether  labeled  or  non- labeled  plots  are  desired. 

13.  Accuracy:  Single  precision 

14.  Caution  to  users : 

Seismogram  must  be  requested  in  the  same  order  that  they 
occur  on  the  input  tape. 

15«  Equipment  configuration:  9  COOP,  23371 , HAWKINS , 
I/3/0/2/4/S/56/57/E/  24=54,  60,5000. 

16.  References:  See  me  and  write-ups  of  Z24  UES  SUBSET  and 
PLOTPACK. 


BLACKY 2 


1102  =  if  1  -  one  channel  per  request 

if  2  -  two  channels  per  request 

NN  =  h  total  number  of  filter  coefficients  to  be  designed. 

0  is  no  filter  is  to  be  designed. 

NDEC  =  Decimator  index,  0  or  1  takes  every  pt,  2  every  other 

point,  etc. 

NDEG  ■  removes  trend  of  degree  NDEG  (usually  zero) 

INTAPE  =  Input  tape  -  always  3 

LAGS  =  Number  of  lags.  LAGS £501 

If  1102  =  1 

ISW1  2  -  semilog  plot  of  power  spectra;  1  -  no  plot 

ISW2  2  -  linear  plot  of  autocorrelation;  1  -  no  plot 

ISW3  2  -  linear  plot  of  power  spectra;  1  -  no  plot 

ISW4  -  always  1 

If  1102  =  2 

ISW1  2  -  semilog  plots  of  power  spectra  of  x  +  y;  1  -  no  plot 

ISW2  2  -  linear  plots  of  autocorrelation  x,  autocorrelation  y, 

cross  correlation  xy,  autocorrelation  yx;  1  -  no 
plot 

ISW3  2  -  linear  plots  of  power  spectra  x,  power  spectra  y, 

co-spectra,  qua-spectra;  1  -  no  plot 

ISW4  2  -  linear  plots  of  coherency  and  phase;  1  -  no  plot 

ICODE  Plot  parameter 

IPTAPE  Plot  tape,  always  2. 

*  if  1  ■  Subset  type  input  tape  with  label 

if  2  -  Subset  type  input  tape  with  label  read  on  some 

preceding  request. 


IAB 


e.g.  -  A  subset  tape  contains  a  label  followed  by  4 
channels  (1)  to  get  power  spectra  only  of 
each  channel 

CARD  1  1102  =  1,  IAB  =  1 


I 


FORMATS 


t 


Card  1  (1613,  3F10 

.1) 

Column 

Data 

• 

3 

1102 

1  or  2 

4-6 

NN 

7-9 

NDEC 

10-12 

NDEG 

i 

13-15 

INTAPE 

-  always 

3 

16-18 

LAGS 

1 

21 

ISW1 

1  or  2 

24 

ISW2 

1  or  2 

27 

ISW3 

1  or  2 

30 

ISW4 

1  or  2 

1 

31-33 

ICODE 

plot  parameter 

j 

1 

36 

IPTAPE 

-  always 

2 

• 

39 

IAB 

1,  2,  or  3 

I 

1 

42 

NSAX 

1  or  2  (usually) 

• 

43-45 

NSKIP 

j 

48 

NFIL 

0  or  1 

49-58 

RANGE 

plot  parameter 

59-68 

CONI 

69-78 

CON  2 

Following 

card  is 

needed  only 

if  NN  0 

FORMAT 

(5F10.3,  110) 

1-10 

FB 

11-20 

FC 

21-30 

FT 

31-40 

SP 

41-50 

RANGE 2 

51-60 

IC2 

Following  card  is  needed  only  if  IAB  *  3 


(16X,  15,  F7.1,  5A8) 


! 


17-21  N 

22-28  DT 

29-68  ID 

COOP  Card 

Input  -  unit  3  -  (subset) 

Output  -  unit  2  -  plot  tape 
unit  4  -  save  tape 

Input  parameter  identical  to  BLACKY 2  with  one  addition 
R  -  used  in  smoothing 


CARD  2  1102  =  1  IAB  =  2 
CARD  3  1102  =1  IAB  =  2 
CARD  4  1102  =  1  IAB  =  2 


(2) 

To  get  cross  correlations  of  channels  1+2  and 

3+4 

f 

1 

CARD  1  1102  =  2  IAB  =  1 

CARD  2  1102  =2  IAB  =  2 

To  get  cross  correlations  of  a  pair  of  channels,  both 
channels  must  be  under  same  subset  label. 

If 

IAB  =  3  -  input  tape  contains  no  label 

NS  AX 

If  NSAX  =  1  power  spectra  is  written  nn  ssuo  t-apa  r.ri  4-v, 
Subset  type  label,  1102  spectra  and  an  and  of  file. 

If  NSAX  =  2  period  and  spectra  only  are  written  on  tape 

NSKIP 

-  number  of  binary  records  to  skip  on  input  tape.  Label 
a  separate  record 

NFIL 

-  If>0  input  data  is  filtered  with  last  filter  setting 
read  in. 

RANGE 

-  plot  parameter  for  linear  plots 

- 

CONI 

-  If  ^  0  each  point  in  first  channel  is  multiplied 
by  CONI 

CON  2 

If  1102  =  2  and  C0N2  ^  0  each  point  in  second  channel 
is  multiplied  by  C0N2 

SEISMIC  DATA  LABORATORY 
Alexandria,  Virginia 
Digital  Computing  Section 


A.  IDENTIFICATION 
Title:  AVGSPEC 

COOP  Identification:  Z133  AVGSPEC 
Category:  General 
Programmer :  S .  Kauder 
Date:  19  April  1967 

B.  PURPOSE 

To  average  spectra  computed  by  program  BLACKY 2 . 

C.  USAGE 

1*  Operational  procedure:  This  is  a  FORTRAN- 6 3  main  program. 
The  binary  input  tape  containing  the  spectra  to  be  averaged  is 
generated  by  BLACKY2.  Any  number  of  jobs  on  the  same  or 
different  tapes  may  be  run  at  one  time.  The  program  stops  on 
an  end-of-f ilfc  on  the  card  reader.  The  plot  tape  (#4)  is 
written  with  program  control. 

2.  Parameters : 

Card  1  FORMAT (4110) 

Col  Description 

I- 10  NS  number  of  spectra  to  be  averaged 

(maximum  of  25.  Maximum  number  of 
points  in  spectrum  is  501) 

II- 20  NLOST  number  of  spectra,  if  any,  to  be 

skipped  on  input  tape. 


If  NN>  0 

FB 

FC 

FT 

SP 

RANGE 2 
IC2 


the  following  paramet  ?rs  are  needed 

beginning  frequency  (cps) 
roll-off  frequency 
cut-off  frequency 

sampling  rate  of  data  (samples/second;,  if  blank 
set  to  20 

plot  parameter  for  filtered  data 

plot  parameter  for  filtered  data.  DATA  is  plotted 
if  IC2  >0. 


IF  IAB  -  3  The  following  cards  are  needed 

N  number  of  points  per  channel 

DT  sampling  rate 

ID  anY  40  character  identification 

Filter  setting  needs  to  be  read  in  only  one  time.  All 
following  requests  with  NFIL  =  1  will  be  filtered  with  this 
setting.  Filter  may  be  changed  on  any  request  by  setting  NN.*  0. 


BLACKY 4 


FORMATS  are  identical  to  BLACKY2  with  following  exception 
Card  1  1613,  4F5.1 


Column 

1-48 

same  as  BLACKY 2 

49-53 

RANGE 

54-58 

CONI 

59-63 

CON  2 

64-68 

R 

BLACKY 4  -  program  control  of  plotting  and  writes  BCD  label  before 
spectra 


Col  Description 

21-30  INTAPE  input  tape  unit.  If  INTAPE  is  zero,  in¬ 

put  unit  is  set  to  3.  (All  input  units 
must  be  indicated  on  COOP  card) 

31-40  NCODE  Determines  length  of  plot  (x)  if  NCODE* 

1,100  pts  per  inch  are  output. 

Second  card  is  needed  only  when  NLOST  is  nonzero. 


Card  2  FORMAT  (25A3) 

1-3  LOST(l)  Channel  ID  of  spectra  to  be  skipped. 

4-6  LOST (2) 

7-9  etc. 

3.  Space  required:  14000  locations. 

4.  Temporary  storage;  none. 

5.  Printouts:  An  appropriate  heading  and  the  average  values 
are  printed. 

6.  Error  returns:  none. 

7.  Error  stop:  none. 

8.  Input  and  output  mountings:  The  input  binary  tape  is  gen¬ 
erated  by  BLACKY 2.  The  input  unit(s)  is  assigned  by  the 
user  by  INTAPE,  (see  section  C-2).  The  output  plot  tape 
is  #4. 

9.  Input  and  output  formats:  See  section  C-2. 

10.  Selective  jump  and  stop  settings:  none. 

11.  Timing: 

12.  Accuracy:  N/A 

13.  Caution  to  user:  none. 

14.  Equipment  configuration: 


15. 


Reference:  Program  BLACKY 2 


MEMORANDUM 


TO: 

FROM: 

SUBJECT: 

DATE: 


Users  of  Avgspec 
E.  F.  Lutz 

Addition  to  Avgspec  Program 
February  15,  1968 


Avgspec  has  been  modified  to  accept  a  variable 
NCODE  to  determine  the  length  of  plots  (X-AXIS) .  The  card 
format  is  now  (4110)  with  columns  31-40  containing  the  value 
of  NCODE. 

An  NCODE  of  1  equc.ls  100  ppi,  2  equals  50  ppi,  etc. 

The  NCODE  value  is  printed  on  the  output. 


EFL: car 


APPLN1'  IX  i: 

Calibration  parameters 


Calibration  for  TF0-J7  on  05  March  1968 


Z01 

13.015  ctr/ni.i 

Z02 

14.959 

Z03 

13. 121 

Z04 

14.403 

Z05 

14. 34 J 

Z06 

14.280 

Z07 

14.090 

Z08 

14.339 

Z09 

14.156 

Z10 

15.569 

Zll 

16.368 

Z12 

15.107 

Z13 

12.937 

ZJ4 

12.194 

Z15 

13.439 

Z16 

15. 490 

Z17 

*14.156 

Z 18 

*14.156 

Z19 

15.906 

Z20 

*14.156 

Z21 

13.  642 

Z22 

10. 304 

Z23 

<13. 110 

Z24 

11.797 

Z25 

*14.156 

Z26 

16.951 

7.21 

16.896 

Z28 

12.463 

Z29 

12.176 

Z30 

12. 195 

Z31 

14.207 

Z  32 

13. 668 

Z33 

13.958 

Z34 

17.880 

Z35 

13. 089 

Z36 

16.734 

Z37 

*14. ] 56 

y  -  i*_ 

“"2,2 
4*  f  hi 

Y  **  Displacement  in  MICRONS 
G  =  Motor  Constant  in 
Newtons/Ampere 
I  “  Current  in  u  amp 
f  *  Frequency  in  CPS 
m  ■  MASS  of  seismometer  in  kg 

For  JM  seismometer  at  leps 
47r2F2m  -  710.6 


*AVG  CAL. -No  CAL  AVAILABLE 
FOR  THESE  SENSORS. 


Him  »«tn«  41  ••  II 

It  III  'l»t  or  •  »l 

—  |»  IT  II. |  IIIHMI  I  IS 

|k*Ut  fill  NO.  If  71«f 


APPENDIX  F 
Sample  output 


•  1  T 
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••N 100  2 

14324 

CM 

ttAWWf 
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S/u«9 
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•  n 

.44 

12.  »4 

44,41 

3.11 

.4* 

•2 

11. *1 
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.32 
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44,14 

2.4* 

.4* 

•  I 

•  •*« 

.44 

11.12 

44,42 

3.44 

,42 
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44,44 
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•» 
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.24 

14.44 

3f.fl 
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.*1 

•4 

Oftitku 

1/ 
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•91 

.25 

14,34 

43.24 

2.',* 

.4* 

•f 

12. or 

•21 

•♦I 

24.34 

41.24 

2.** 

,2« 

If 

lt»»2 

•23 

.*4 

12.41 

4|.«l 

2. *7 

.4* 

1# 

-  _  io*rf 

•9* 

•  *2 

12.4# 

41.21 

-  3.43— 

11 

«•<« 

•  1« 

f./J 

44.99 

4.11 

.4* 

12 

If. 44 

•  9* 

.•4 

12.47 

43.24 

3.12 

lot* 

13 

f  .42 

•  94 

.»f 

13.14 

44.11 

2.44 

1*11 

14 

T.23 

.i* 

.24 

12.41 

34.44 

3.4* 

loll 

l* 

4. *7 

•i» 

.«• 

11.  J2 

34.21 

2.4* 

.»! 

14 

*•19 

•  if 

•  99 

If  4  *4 

Jl  all 

-  2*9*4 

,94 

IT 

4.  >i 

•  1* 

•  9 1 

11.22 

14.44 

2.43 

.41 

11 

f  .>7 

•  92 

.*• 

12.44 

St  •  24 

2*lt 

.2* 

If 

14.22 

•  21 

.•1 

12.44 

33. «f 

2  •*• 

•  it 

24 

If  .4f 

•  37 

1.12 

14.44 

14. 11 

3.43 

!  * 

21 

11.34 

•  9? 

l.»T 

4.24 

42.31 

T.Tl 

.22 

22 

.  .  lllHtii 

22 

11. If 

•  9t 

1.44 

14.41 

40.31 

3.*3 

.41 

24 

11.34 

.22 

1.34 

2.12 

42.11 

4.T4 

.•« 

2» 

•  .»! 

•  IT 

.*2 

11. *1 

41. 44 

1.44 

l,|* 

24 

«.4» 

.  1  T 

.«f 

14.21 

14,31 

?••* 

t*2f 

2T 

4.13 

•  22 

.31 

11.24 

42. ft 

3.44 

1.31 

_ at 

- i,a 

«?? 

,14  ... 

2f 

4.34 

•  1* 

>7 

4.21 

14,44 

2.f  * 

— - 1  T*+ - 

|,3* 

34 

1.24 

.11 

.3* 

1.1> 

11*42 

3,41 

|,l/ 

31 

4.34 

.24 

.3* 

4,44 

31.42 

1.4* 

»• 

32 

4,«7 

•  9* 

./ 4 

4.44 

40.34 

2. *• 

:*r 

32 

4.32 

•if 

.*1 

10.42 

34.44 

3.24 

.44 

J4-. 

_ tt'l 

,t». 

.»» 

14  t  9* 

40.14 

•  fl 

3* 

HWl 

•  9* 

.*• 

14. f4 

42.  »3 

2.4* 

,tt 

31 

14. 'I 

.3* 

1.24 

.  12.31 

41, *4 

1.41 

.1* 

IT 

OUifftl 

IUN 

l|2«2» 

T.Tl 

2*.l* 

«2|.*2 

1244.42 

11*. 3l 

Mf*N 

••I* 

.23 

12.3* 

»  2  ,  T* 

3.3* 

•  •• 

f  .24 

•ft 

ill 

13. 2» 

1*4.2* 

.*» 

•to 

3.44 

•  ft 

.«• 

3.44 

43.T* 

.*• 

- 

U*»UN 

Ilf} 

•  4» 

•  If 

2.43 

22,1* 

3. if 

•••UN 

4,4| 

•  ft 

.!» 

44,42 

1*4,23 

3,41 

»U*fUN 

1.2* 

•  «» 

.1* 

f.4l 

24.2# 

i.*r 

••••UN 

4.»f 

•  fit 

.14 

41.41 

1*4.1* 

3.40 

•IONIC  IU(| 
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• 

NO  Ilk  I0SI 

4/N  1041 
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APPENDIX  G 

JM  frequency  response 


SHORT  PERIOD  -  (JOHNSON  MATHESON) 

Seismometer  free  periods  1.25  sec 
Seismometer  damping:  0.510 

Galvanometer  free  period:  0.333  sec 
Galvanometer  damping:  0.650 


T 

f 

« 

Norm 

Max 

Norma  1  i  Min 

5.0 

0.2 

10 

0.0120 

0.0132 

iTuiiiiaii  icu  nin 

0.0108 

2.5 

0.4 

7.5 

0. 0950 

0.102 

0.0878 

1.25 

0.8 

5 

0.685 

0.720 

0.651 

1.0 

1.0 

0 

1.00 

1.00 

1.00 

0.67 

1.5 

5 

1.52 

1.60 

1.44 

0.5 

2.0 

5 

1.90 

2.00 

1.  80 

0.33 

3.0 

7.5 

2.12 

2.28 

1.96 

0.25 

4.0 

12 

1.87 

2.09 

1.64 

0.167 

6.  0 

20 

1.15 

1.38 

0.920 

0.125 

8.  0 

25 

0.64 

0.  800 

0.  480 

0.1 

10.0 

25 

0.370 

0.463 

0.278 
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-  Beamforming  of  the  TFO-37  array  reduces  the  rms  of  the  noise 
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